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In this paper we present a fully ab initio Hartree-Fock 
approach aimed at calculating the static structure factor of 
crystalline insulators at arbitrary values of momentum trans- 
fer. In particular, we outline the computation of the incoher- 
ent scattering function, the component of the structure fac- 
tor which governs the incoherent x-ray scattering from solids. 
The presented theory is applied to crystalline LiF to obtain 
benchmark Hartree-Fock values for its incoherent scattering 
function. Benchmark theoretical values such as this, can be 
combined with the experimentally measured static structure 
factor, to understand the influence of electron correlation ef- 
fects on cohesive properties of solids. 



I. INTRODUCTION 

In order to obtain an ab initio understanding of the 
electronic structure of solids, it is essential tOj-nnderstand 
the nature of electron correlations in themtl. Correla- 
tion being a real-space phenomenon, one possible way of 
visualising it is through the density-density correlation 
function defined as 



S(r',r) 



1 



($|p(r')p(r)|$) 



(1) 



where $} denotes the many-particle wave function of the 
system, N is the total number of electrons in the system 
and p(r) — Yli^i — r «) i s the density operator with r, 
being the coordinates of the i-th electron. One can easily 
show 



S(r',r) =S(r' -r) + (N -l)g(r',r) 



(2) 



where <?(r', r) is the electron pair-correlation function de- 
fined as 

5(r '' r) = ^ 6{T ' ri)6{r ~ Fj)m • (3) 

Above p(r) denotes the electronic charge density. Pair- 
correlation function represents the probability that when 
one electron is observed say at point r', another electron 
will he found in a characteristic volume 1/(Nq — 1) located 
at ra. Clearly it can be used to quantify the so-called 
"exchange-correlation holeli associated with an electron 
in a many-electron systerna. If we compute the Fourier 
transform of the density-density correlator, we obtain the 
static structure factor 5(Q) defined by 



5(Q) = ^ J drdr'e iQ < T - r,) S{v',v) 



(4) 



If we assume that the system under consideration is a 
crystalline system with N unit cells, each of which has 
Z electrons in it so that Nq = NZ, we can decompose 
5(Q) as a sum of a "coherent" and an "incoherent" part 

S(Q) = ^<5 Q , G |F(Q)| 2 + S inc (Q) , (5) 

where G is a vector of the reciprocal lattice and the form 
factor F(Q) defined as 



No 



F(Q) =^($|e lQ - r ^|$) 



(6) 



can be easily seen to be the Fourier transform of the 
charge density of the system, while 

S inc (Q) = — ( £ ($|e^-^-^|$) - |F(Q)| 2 ) (7) 



j,k=l 



is referred to as the incoherent scattering function in the 
literature. It is intuitively obvious that being the ex- 
pectation value of a two-electron operator, S'm C (Q) will 
be sensitive to electron correlations in the crystal, while 
F(Q), which is a one-electron operator, should be rela- 
tively insensitive to such effects. It is easy to verify that 
•Sinc(Q) satisfies limiting conditions 



Hm S inc (Q) 







and 



lim 5m C (Q) = 1 



(8) 



(9) 



One can perform the measurement of the static struc- 
ture factor of a many-electron system in a variety of ex- 
periments such as electron scattering and x-ray scatter^, 
ing. For isplated atoms and molecules both electronscl 
and x-rayaj are frequently used for such measurements, 
however, for crystalline sytema., x-ray scattering appears 
to be the method of choiceQH. In such measurements, 
the quantity Q in equations above is identifed with the 
momentum transferred by the incident particle (electron 
or photon), to the many-electron system under investi- 
gation. Keeping in mind the relationship between the 
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static structure factor and the pair-correlation function 
(cf. Eqs. (||) and ([IJ) ), it is thus possible to obtain the 
pair-correlation function from these measurements. How- 
ever, in what follows, we will devote exclusively on the 
x-ray scattering based experiments. The coherent x-ray 
scattering (Bragg scattering), i.e. when the momentum 
transfer Q is equal to one of the vectors G of the re- 
ciprocal lattice, is governed predominantly by the form 
factor F(Q). However, by concentrating on the mea- 
surements corresponding to those values of momentum 
transfer which are not equal to any reciprocal lattice vec- 
tor, one can — according to Eq.(^|) — directly measure the 
incoherent scattering function. The experiments which 
concentrate on this region of Q correspond to incoher- 
ent x-ray scattering. For the case of incoherent scatter- 
ing of x-rays from a crystalline solid at finite tempera- 
tures, assuming that the energy of the incoming x-rays 
is much higher compared to the binding energies of the 
constituent electrons, but still low compared to the rest 
energy of the electron moc 2 (so that the relativistic ef- 
fects can be neglected), the scattering cross section, for 
a pp)^d composed of light elements, can be approximated 
as 



3n * rlN 



■f- ) (e~ • e) 2 Si, 

«0 ' 



:(Q) + S T DS(Q) 



, (io) 



where ro is the classical electron radius, ko and k are the 
wave vectors of incoming and outgoing photons, e~o and 
e are the corresponding polarization vectors, Q = k — k 
is the momentum transfer while Stds(Q) is the struc- 
ture factor due to the thermal diffuse scattering (TDS) 
caused both by the thermal, and the zero-point vibra- 
tions of the lattice. Thus Stds{Q) quantifies the con- 
tribution of phonons to the x-ray scattering, and can be 
computed by-taking one-phonon, and higher order terms 
into accountOu. Therefore, by measuring the incoherent 
x-ray scattering cross section for different values of the 
momentum transfer Q, combinded with the knowledge of 
Stds(Q), one can, using Eq. (|io|), extract the incoher- 
ent scattering scattering function Sjnc(Q) of the system 
under consideration. nn 

Sacchetti and coworkera3u have been the proponents 
of using the incoherent x-ray scattering to measure the 
static structure factor of crystalline compounds, for its 
subsequent use in the analysis of electron correlation ef- 
fects. They have performed a series of accurate mea- 
surements rpf the static sturcture factor of the .-metallic 
system BeQ and the covalent system diamondQ to ob- 
tain their pair-correlation functions, and analyzed vari- 
ous contributions to the ground state energies of these 
compounds. In their latest experiment performed on 
crystalline LiF, they have, for the first time, subjected 
an ionic system to a similar analysis!! However, in order 
to quantify the contribution of electron correlation effects 
to the experimentally measured static structure factor in 
such experiments, benchmark Hartree-Fock (HF) results 
for the quantity are needed. It is the purpose of this 



paper to present a formalism using which one can per- 
form such benchmark HF calculations within an ab initio 
framework. Indeed, Calzuola et alB, by comparing their 
experimentally measured values of the static structure 
factor, to the benchmark HF values presented here, have 
estimated the correlation contribution to the cohesive en- 
ergy of LiF. Therefore, the aim of these calculations is not 
to explain the experimental data, but rather to provide a 
theoretical reference, with respect to which the correla- 
tion effects can be quantified in the experimentally mea- 
sured quantities. The formalism for computing the static 
structure factor presented here is based on a Wannier- 
functionJpascd ab initio HF approach developed recently 
by ugilJ'tJ. The approach has since been applied to fGip* 
pute the ground state pro-parties of a number of ionidi3'E3 
and covalent cpmpounds£-lt3, including the form factor 
(F(Q)) of LiFlII 

Remainder of the paper is organized as follows. In sec- 
tion |n] we describe our formalism for the ab initio evalu- 
ation of the incoherent scattering function within an HF 
approach. An explicit formula is presented which rep- 
resents S l j nc (Q) in terms of the Wannier functions of an 
infinite crystal. Our numerical results for LiF are pre- 
sented in section III which arej-jcompared to the experi- 



mental results of Calzuola et alB for the same compound 
Finally our conclusions are presented in section IV. 



II. THEORY 

Here we outline the evaluation of the incoherent scat- 
tering function «Si„c(Q) for an infinite crystalline insu- 
lator within an ab initio restricted Hartree-Fock (RHF) 
approach. Although, we are not aware of such a prior 
calculation for an infinite solid, we note that ab initio 
calculators are performecLon-a routine basis on isolated 
atomsB~E3 and moleculesEm, both at the HF, and the 
correlated level. However, for condensed-matter systems, 
perhaps because of practical difficulties associated with 
providing a wave-function-based ab initio description of 
an infinite system, such calculations are either performed 
assuming a jellium model for the electrons of the sys- 
tcmou, or within the framework of the density functional 
theory which often involves phenomenological approxi- 
mationscJ. 

We assume that the compound under consideration is 
a closed-shell crystalline system whose RHF ground state 
can be described by n c doubly occupied Wannier func- 
tions per unit cell, so that Z — 2n c . If we use Greek 
indices a, j3 etc. to denote the Wannier functions lo- 
calized in a given unit cell, the RHF wave function of 
the infinite crystal can be described as a Slater determi- 
nant composed of the infinitely many Wannier functions 
{|a(Rj));a = l,n c ;j = 1,N}, where |a(Rj)) denotes 
the a-th Wannier function of a unit cell located at the 
position given by the vector Rj of the lattice. The afore- 
mentioned Wannier functions are assumed to form an 
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orthonormal set 

(a{R i )\p(Il j )) = 5 a p6 ij , (11) 

and Wannier functions localized in different unit cells are 
translated copies of each other 

\a(Ri + Rj)) = T(R l )|a(R J )), (12) 

where the operator 7~(Rj) represents a translation by lat- 
tice vector Ri. The theory and several applications of our 
approach, which directly obtains the RHF Wannier func- 
tions of a crystalline insulator feriie or covalent), have 
been presented in several paperst^L3. 

If we use the standard formula for the expectation 
value of a two-particle operator with respect to a Slater 
determinant^, one can, after some algebraic manipu- 
lations, show that in the Wannier representation the 
RHF expression for the incoherent scattering function 
(cf. Eq.(0)) is given by 

„ n c N 

S inc (Q) = l-- Yl El^)l elQ >(°)>l 2 ' (13) 

a,P=l i=l 

where \a(0)} represents a Wannier function localized in 
the reference unit cell. Eq. (|l^) constitutes the key for- 
mula of this work, a detailed derivation of which is pre- 
sented in the appendix. Since the Wannier functions in 
our computer code are represented in terms of Gaussian 
lobe-type localized basis functions, it is possible to write 
down analytic expressions for the matrix elements needed 
to evaluate Sj nc (Q) according to the expression above. 
The lattice sum over lattice vectors R^ involved in Eq. 
( |l3| ) decreases rapidly as one moves away from the refer- 
ence cell, and is terminated once the convergence within 
a given threshold is achieved. 

The restriction of the present approach to insulators 
stems from our use of Wannier functions as the single- 
particle orbitals, rather than the conventional Bloch or- 
bitals. It is easy to see that the same theory can be easily 
extended to metallic systems if one were to express the 
many-body wave function of the solid in terms of Bloch 
orbitals. In that case, of course, the real-space sum (cf. 
Eq.(|l3|)) will have to be replaced by an integration over 
the Brillouin zone. We will present this generalization in 
a future paper. This will be particularly useful in light 
of the future experiments which Sacchetti and coworkers 
are planning on metallic systemsLl 

III. RESULTS AND DISCUSSION 

In this section we present the results of our HF calcu- 
lations of the incoherent scattering function at different 
values of the momentum transfer, and compare our re- 
sults to the experimental ones. It is intuitively obvious, 
however, that HF structure factors can only describe the 
experimental results qualitatively — in order to obtain a 
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FIG. 1. ZSi„c(Q) plotted as a function of the momen- 
tum transfer Q in the direction [100]. Solid line represents 
the HF theoretical results of this work, while the dashed line 
represents the experimental results of Calzuola et al.u. The 
experimental data was confined to the values of Q ranging 
from 0.517 a.u. to 7.622 a.u. 

better quantitative description theoretically, inclusion of 
electron correlation effects is essential. Nevertheless, in 
our opinion, the comparison with experiments is very in- 
structive, because one can, in a rather pictorial way, see 
the successes and failures of the HF approximation in 
describing the physics of weakly-correlated systems. 

The basis set used to represent the Wannier functions 
in our calculations was the lobe representation of the ba- 
sis set proposed by Prencipe et al. in their Bloch-orbital- 
based HF study of the structural properties of LiFE3. 
The basis set consisted of contracted Cartesian Gaussian- 
type basis functions and was of [4s, 3p] type for the fluo- 
rine atom, and [2s, lp] for the lithium atom. For further 
details pertaining to the exponents andj-the contraction 
coefficients we refer to the original workEi Details deal- 
ing with the lobe representation of the Cartesian Gaus- 
sian basis functions can be found, e.g, in our previous 
papertil. We also examined the basis-set dependence of 
our results on Si nc (Q) by performing calculations with 
larger basis sets which also included d-type basis func- 
tions on F atom, however, we did not observe any sig- 
nificant change in the results. Thus we believe that our 
results on Sj nc (Q) presented below are fairly accurate. 

In the theoretical calculations the observed face- 
centered cubic (fee) structure was assumed for the com- 
pound. The reference unit cell was taken to be the prim- 
itive cell with the F atom at (0, 0, 0) position and the Li 
atom at (0,0, a/2), where a is the lattice constant. For 
the lattice constant, the room temperature value of 4.02A 
was used. 

The results of our calculations are are summarized in 
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table | and Fig. [l] which present ZSi nc (Q) as a func- 
tion of the momentum transfer Q, where Z — 12 for 
the case of LiF. Direction of the momentum transfer for 
both the theory and the experiment was along [100] direc- 
tion. Thcoreitcal HF values of Si nc (Q) were calculated 
with the magnitude of the momentum transfer Q ranging 
from 0.0 to 15.0 atomic units (a.u.). For the experimental 
data, Q ranged from 0.517 a.u. to 7.622 a.u. It is quite 
clear from Fig. |l| that HF theory is in good qualitative 
agreement with the experimental result, which is a man- 
ifestation of the fact that LiF is a weakly correlated sys- 
tem. In order to quantify the correlation effects, we define 

the quantity E(Q) = S — ( JE~^p q) x 100, which clearly 

measures the percentile contribution of electron corre- 
lation effects to the experimentally measured Sj nc (Q), 
using the HF values presented here as the benchmark 
reference. E(Q) is plotted in Fig. ||, as a function of 
Q. For the smallest value of the momentum transferred 
measured Q = 0.517 a.u., the correlation contribution 
is 88.2%. With the increasing momentum transfer the 
correlation contribution decreases rapidly staying in the 
range 10.0 % — 20.0 % from Q = 0.724 a.u. to Q = 2.059 
a.u.. From Q = 3.664 a.u. onwards the uppper bound for 
the correlation contribution is approximately three per- 
cent, while most of the points are in one to two percent 
range. We also see some oscillations in the E(Q) as a 
function of the momentum transfer, which may be due 
to the experimental uncertainties. However, the general 
trend in E(Q) as a function of the momentum transfer 
Q is clear — the contribution due to the correlation ef- 
fects decreases with the increasing momentum transfer. 
This trend is also observed in the calculations involving 
free atoms and molecules where the HF calculations for 
small values of momentum transfer always Oiier-estimate 
Sinc(Q) as compared to the correlated onestSTo. This 
trend can be understood as follows. HF theory, because 
of a lack of correlations in it, will always overestimate 
the pair-correlation function g(r',r) and consequently 
5i nc (Q). Since for small values of momentum transfer 
Q, one can only probe the valence electrons, the main 
contribution to SV, lc (Q) will also naturally come from 
these electrons. However, it is the valence electrons for 
which the correlation effects are quantitatively the most 
important, and their neglect in the HF approach leads 
to relatively large deviations, as compared to the experi- 
mental values, for smaller values of momentum transfer. 

In our calculations we did not observe any significant 
anisotropy in 5j nc (Q) with respect to the direction of Q. 
We performed the same set of calculations for Sj nc (Q) 
for momentum transfer directions [110] and [111] as well, 
however, the difference in the results compared to [100] 
direction was always less than 1.0 x 10 -5 . This result can 
also be understood on intuitive grounds as the charge 
density in LiF is fairly isotropic, therefore, one would 
not expect the incoherent scattering function to show 
any anisotropy. Finally it is clear from both table || and 
Fig. [I] that our HF results on 5i„ c (Q) approach the cor- 



TABLE I. Comparison of the Hartree-Fock incoherent 
scattering function (Z Smc{Q)) computed in this wqck with 
those measured in the experiment of Calzuola et alfl at se- 
lected values of momentum transfer Q. The momentum 
transfer was along the [100] direction in both the experiment 
and the theory. 
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FIG. 2. Relative correlation contribution E(Q) to the mea- 
sured Sirac(Q), plotted as a function of the momentum trans- 
fer Q. The momentum transfer direction was [100]. See text 
for details. 



larging the many-particle ground-state wave function by 
considering virtual excitations from the space of the oa- 
cupied Wannier functions to that of the virtual onesEa. 
The approach was demonstrated by computing the cor- 
relation cuantributions to the total energy per unit cell of 
bulk LiHE£l. However, the generalization of the approach 
to compute the correlated expectation value of an oper- 
ator other than the Hamiltonian is far from trivial, and 
will be the subject of a future investigation. 
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rect limiting values in both the low momentum-transfer 
region (cf. Eq.(|l)) and the high momentum-transfer re- 
gion (cf. Eq.ffif- 



IV. CONCLUSIONS 

In conclusion a wave-function-based fully ab initio ap- 
proach has been presented, using which, one can compute 
the static structure factor of a crystalline compound at 
arbitrary values of the momentum transfer, at the HF 
level. The formalism was applied to the case of crys- 
talline LiF, and benchmark values were obtained for its 
incoherent scattering function. These values were used 
in the analyssis of a recently performed incoherent x-ray 
scattering experiment on LiF to quantify the electron cor- 
relation effects, in general, and to predict the correlation 
contribution to its cohesive energy, in particular!! In case 
such experimental measurements are performed on other 
insulating compounds, one can use the formalism pre- 
sented here to perform similar benchmark calculations 
on those systems as well. The present version of the 
theory is restricted to insulating systems because of its 
Wannicr-function-based formulation, however, in light of 
the planned future experiments on metallic systemstj, we 
do intend to develop an ab initio HF formalism meant for 
computing Si, lc (Q) for gapless systems, as well. 

Although the aim of the present HF formalism was not 
to explain the experimental data, but rather to facilitate 
its analysis, it is still of interest to include electron cor- 
relation effects theoretically. Indeed, recently, we have 
generalized our Wannier-function-based approach to in- 
clude electron correlation effects by systematically en- 



APPENDIX: 



Our aim in the present section is to present a derivation 
of Eq.([t3]) of the text which is an RHF level expression 
for S'm C (Q) in terms of corresponding Wannier functions. 
Eq.(0) which defines Si nc (Q) involves the expection value 
of the operator 



A(ri,r 2 , . . .,r No ;Q) 



N 

= eiQ ' [ri 



(Al) 



Although operator X(ri, ra, . . . , rjv ; Q) is a sum of two- 
electron terms, however, unlike other similar operators 
such as the Coulomb interaction operator, the sum in 
Eq.(Al) does contain the term where j = k. Therefore, 



in order to utilize the well-established formulas for ma- 
trix elements-, of two-electron ope rators between Slater 
determinants, we rewrite Eq.(Al) as 



^(ri,r 2 , 



,rAT ;Q) =iV +y(ri,r 2 ,...,rjv o ;Q) 

+y*(r 1 ,r 2 ,...,r Wo ;Q) , (A2) 



where the first term on the right-hand side corresponds 
to j = k terms of the sum in Eq.(Al) and the operator 
Y is defined as 



y(ri,r 2 ,...,r iVo ;Q) 



iQ-(r,-r fc ) 



(A3) 



Y* , which represents the complex c onj ugate of operator 
Y, can be easily deduced from Eq.(|A3|). It is clear that 
Y (and hence, Y*) as defined above are traditional two- 
electron operators. Utilizing the well-known formula for 
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the expectation value of a general two-electron operator 
G = J2j <k g( r j, ttJ with respect to a single Slater deter- 
minant state |\I/)e3 

(*|G|*) = i5>%((r 1 ,r 2 )|ab) 



a,b 



l^2( ba \g((*i, r 2)\ab) , 



(A4) 



where a), \b) are the orbitals constituting the Slater de- 
terminant, one can easily get 



N 

^ E 

a, 0=1 Ri,R 3 =l 



{2(a(R, 



a(Ri)> 



x^Ole-^'l^)) " (/3(R,)| e ,;Q ' r |«(R 4 )) 



x(a(R0|e- lQ - r |/3(R,)} , 



(A5) 



where |$) represents the single-Slater-determinant RHF 
ground state of the crystalline insulator under consider- 
ation. As discussed in Sec. |l[ it is expressed in terms of 
Wannier functions {\a(Rj));a = l,n c ;j = 1,N}, where 
2n c is the number of electrons per unit cell and N is 
the total number unit cells considered. In order to ar- 



rive at the right-hand side of Eq.(A5) spin summations 
have been performed, and the first term represents the so- 
called "direct" contribution, while the second term rep- 
resents the "exchange" contribution. It is clear that both 
the direct, as well as the exchange terms are expressed as 
products of two one-electron matrix elements which are 
complex conjugates of each other. Moreover, using the 
defining Eq.(||), one can easily see that, for the RHF state 
considered here, the direct term of Eq. (|A5|) is nothing 
but the product of the form factor F(Q) and its complex 
conjugate, leading to 

(*|y|*> = \\f(q)\ 2 



E 



N 

E 



K^R^Ie^XR,)}! 2 . (A6) 

a, 0=1 Ri,Rj=l 

Finally, using the translational invariance property of 
the one-electron matrix elements (/3(Rj)|e l( ^' r |Q;(Ri)} = 
(/3(Rj — Ri)|e 4< 3' r |a(0)}, and by r earr anging the sum over 
lattice vectors R^ and Rj in Eq.(A6), one gets in the in- 
finite solid limit (N — > oo) 



mm = -if(q)| 2 



E 

a,/3=l R«=l 



^E 



(/Wle^XO))! 5 



(A7) 



If we combine the results of Eqs. (|Al|) , QA2|) , (|A7|), 
and substitute them in Eq.(^), we immediately obtain 
Eq. (|l3|) , valid for a single Slater determinant RHF wave 
function 1$). 
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